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Summary 

A  study  was  made  of  ihe  response  of  a  rotor-bearing 
model  of  a  turbine  engine  rotor  for  two  different  blade 
tip-seal  interference  rub  models.  A  direct-integration, 
transient -response  rotor  dynamics  computer  code  was 
used  for  each  simulation.  The  first  model,  an  abradable 
seal  rub  model,  is  based  on  an  energy-loss-per-unit- 
volume  theory.  This  model  is  applicable  to  blades  coated 
wit  h  an  abrasive  material  rubbing  on  a  ceramic -coated 
seal.  The  second,  a  smearing  model,  is  based  on  a 
Newtonian,  low-Reynolds-number,  viscous  hydro- 
dynamic  theory.  This  model  is  applicable  to  uncoated 
blades  rubbing  on  a  metallic  seal,  which  produces  a  thin, 
molten  metal  layer  between  the  blade  tip  and  the  seal 
substrate. 

The  rotor-bearing  system  used  with  each  rub  model 
consisted  of  a  shafi  with  three  disks  mounted  on  two  stiff 
bearings  that  were  mounted  in  a  squeeze-film  damper 
with  centering  springs.  Two  of  the  disks  were  overhung 
and  the  third  was  centered  with  respect  to  the  bearing 
supports.  The  rotor-bearing  system  had  three  lateral 
bending  critical  speeds  of  130,  150,  and  190  Hz  (7600, 
9200.  and  II  200  rpm). 

The  rotor  system  was  assumed  to  be  balanced  prior  to 
blade  loss  and  operating  at  160  Hz  (9550  rpm).  The  blade 
loss  was  simulated  by  an  instantaneous  application  of  a 
130-^m  (5-mil)  mass  eccentricity  in  only  one  of  the 
overhung  disks.  Each  disk  was  surrounded  by  a  sea! 
shroud  that  had  a  50-^m  (2-mil)  radial  clearance  and  a 
20-MN/m  (100  000-lb/in)  radial  stiffness.  Each  disk  had 
rectangular  blades  with  a  tip  velocity  of  100  m/s  (3800 
in/s),  and  only  one  blade  per  disk  was  rubbing  against  the 
seal  substrate  at  a  time. 

The  rotor  response  was  determined  for  each  rub 
model.  Two  rub  cases  were  simulated  for  each  rub  model. 
In  the  abradable  model  the  energy  loss  parameter  was 
varied.  The  u suits  for  a  minimum  and  maximum  energy 
loss  are  presented  in  this  report.  In  the  smearing  model 
the  viscosity  parameter  was  varied  and  the  results  for  low 
and  high  viscosity  are  also  presented. 

The  rotor  response  for  low  energy  per  unit  volume  in 
the  abradable  model  was  similar  to  the  response  for  low 
viscosity  in  the  smearing  model.  In  both  cases  the  rotor 
swept  out  an  orbital  path  that  stabilized  in  the  forward 
whirl  direction  to  the  same  light  interference  of  5  (0.2 

mil).  Prior  to  steady-state  response  the  amplitude  of  the 
blade  bouncing  against  the  elastically  supported  seat 
substrate  was  small  and  had  a  high  frequency  of  550  Hz 
(33  000  cpm).  The  blade  tip-seal  rub  during  this  time 
interval  was  almost  continuous. 

The  rotor  response  for  high  energy  pet  unit  volume  in 
the  abradable  model  was  similar  to  the  response  for  high 
viscosity  in  the  smearing  model  up  to  the  time  of 
backward  whirl.  During  this  time  interval  the  amplitude 


of  blade  bouncing  against  the  seal  was  larger  and  less 
continuous  than  when  the  corresponding  low  parameters 
were  used  for  these  models.  The  frequency  of  bouncing 
for  each  was  about  400  Hz  (24  000  cpm).  After  25  ms  the 
rotor  went  through  a  transition  from  forward  to 
backward  whirl  for  the  high-energy  and  high-viscosity 
cases.  For  the  high-viscosity  case  and  the  smearing 
model,  the  amplitude  of  backward  whirl  indicated  a  finite 
limit  to  the  orbital  path  of  the  rotor.  For  the  high-energy 
case  and  the  abradable  model  the  amplitude  of  backward 
whirl  grew  without  limit,  similar  to  the  results  of  a  dry 
friction  model. 


Introduction 

In  a  typical  aircraft  gas  turbine  there  are  many 
instances  in  which  rotor  rubs  occur.  Two  of  the  most 
common  are  blade  tip  and  seal  rubs,  which  are  caused  by 
thermal  mismatch,  rotor  imbalance,  high  "g”  maneuver 
loads,  aerodynamic  forces,  etc.  Current  interest  in  fuel 
efficiency  drives  the  engine  design  toward  closer  opera¬ 
ting  clearances  and  thus  increases  the  probablity  of  rotor 
rubs. 

it  is  known  that  rotor  rubs  can  have  an  important 
effect  on  the  rotor  dynamics.  When  a  rotor  rubs  on  the 
case,  a  frictional  force  is  generated  that  can  drive  a  rotor 
to  whirl  in  a  direction  opposite  to  the  direction  of 
rotation  (backward  whirl).  This  frictional  force  is 
relatively  constant  up  to  the  backward  whirl  speed  at 
which  the  rotor  starts  to  roll  around  the  case.  Since  this 
rolling  contact  speed  is  proportional  to  the  rotational 
speed  of  the  rotor  times  the  ratio  of  the  rotor  diameter  to 
the  rotor  clearance,  the  whirl  speed  can  be  hundreds  of 
times  the  rotational  speed  of  the  rotor  and  thus  can  be 
very  dangerous. 

Studies  of  the  interaction  of  a  rotor  with  its  case  (rotor 
rubs)  have  been  reported  in  references  1  to  5.  Reference  1 
concerns  a  steady-state  interaction  between  a  rotor  and  a 
rigid  case  with  friction  at  the  interface  neglected;  and 
reference  2,  a  steady-state  interaction  between  a  linear, 
flexible  rotor  and  case  with  friction  at  the  interface 
included.  References  1  and  2  did  not  consider  the  critical 
transient  situation  in  which  the  rotor  bounces  off  the 
case. 

References  3  to  5  studied  the  transient  situation  in 
which  a  blade  loss  induced  a  rotor  rub.  Reference  3 
neglected  the  frictional  force  of  the  rub.  References  4  and 
5  both  analyzed  the  case  of  a  dry  friction  rub,  that  is,  the 
frictional  force  proportional  to  the  normal  force.  The 
difference  between  references  4  and  5  is  that  r-ferener  4 
used  the  component  mode  synthesis  method  and  refer¬ 
ence  5  used  a  direct  integration  method. 

When  a  rotor  rubs  the  case,  a  complex  thermo¬ 
mechanical  phenomenon  occurs.  Two-dimensional  finite 


element  codes  have  been  written  to  predict  this 
interaction  (ref.  6).  It  is  impractical  to  try  to  incorporate 
this  kind  of  code  into  a  transient  rotor  dynamics  code. 
Therefore  a  model  of  the  rub  phenomenon  is  needed,  it  is 
unlikely  that  the  dry  friction  rub  model  used  in  references 
4  and  5  would  be  adequate. 

The  gas  path  seal  materials  seem  to  be  in  two  classes. 
The  blade-seal  material  combination  of  metal  to  metal  is 
an  older  technology.  A  ceramic  seal  with  blades  coated 
with  an  abrasive  material  is  the  newer  technology. 
Experiments  on  both  combinations  performed  in  steady 
in-house  rub  testers  imply  two  different  rub  models.  The 
metal  blade  rubbing  on  a  metal  seal  casing  leads  to  the 
blade  tip  melting  off  and  smearing  on  the  colder  metal 
seal.  The  abrasive-tipped  blade  rubbing  on  a  ceramic  seal 
leads  to  the  blade  grinding  or  abrading  the  ceramic  seal. 

It  is  the  object  of  this  report  to  incorporate  the 
smearing  and  the  abrading  rub  models  into  the  analysis 
used  in  reference  5  and  to  compare  the  rotor  dynamics  of 
these  two  models  with  that  of  the  dry  friction  model. 


Symbols 

a  blade  tip  thickness 

b  blade  tip  width 

C  radial  clearance  between  blade  tip  and  seal 
Z-!  normal  blade  tip  force 

F 2  tangential  blade  tip  force 

h  thickness  of  molten  metal 

K  radial  stiffness  of  seal 

/  rub  length  on  abradable  surface 

p  pressure  across  liquid  film 

t  time 

U  energy  per  volume  of  material  removed 
u  tangential  blade  tip  velocity 

v  radial  blade  tip  incursion  velocity 
r  radial  displacement  of  blade  tip 

x.y.z  coordinate  variables 

n  viscosity  of  molten  metal 

r  shear  stress  of  molten  metal 


Analysis 

The  rotor-bearing  system  is  assumed  to  be  rotating  at  a 
constant  angular  velocity  with  a  mass  eccentricity  in  one 
of  the  disks.  When  the  radial  displacement  r  of  a  blade 
exceeds  the  radial  clearance  C  between  the  bladed-disi. 
and  the  seal,  normal  F{  and  tangential  Fi  blade  tip  forces 
are  generated  during  rub.  During  a  rub  the  tangential 
blade  tip  velocity  u  is  assumed  to  remain  constant  while 


the  radial  tip  incursion  velocity  v  varies  with  time.  The 
seal  substrate  for  each  rub  mode!  is  assumed  to  be 
elastically  supported  only  in  the  radial  direction  with  a 
stiffness  K. 


Smearing  Rub  Model 

When  a  metal  blade  rubs  on  a  metal  seal  casing,  the 
blade  tip  melts  and  deposits  a  smeared  layer  on  the  colder 
metal  surface  as  the  blade  passes  the  point  of  rub.  At  the 
instant  of  rubbing  a  layer  of  mo’  n  metal  between  the 
blade  tip  and  the  seal  substrate  velops  as  shown  in 
figure  1.  A  tangential  blade  tip  force  is  developed  due  to 
the  blade  moving  with  a  velocity  u  in  a  viscous  medium. 
In  addition,  a  normal  force  is  developed  due  to  the 
relative  radial  velocities  of  the  blade  tip  and  the  seal 
substrate.  These  forces  are  analogous  to  the  forces  in  a 
squeeze-film  bearing.  Assuming  the  ratio  of  blade  width 
to  thickness  b/a  to  be  large,  the  Reynolds  equation  for  an 
infinitely  wide  blade  moving  through  a  liquid  at  a 
constant  velocity  u  and  at  constant  density  and  viscosity 
is  presented  as 


h 3  d*p  _  dn 
12m  3*2  ~  3/ 


(1) 


where  p(x,t)  is  the  molten  metal  pressure  and  h  is  the 
thickness  of  molten  metal  between  the  blade  tip  and  seal 
substrate. 

If  the  boundary  conditions  on  p(x,t)  are 
p(0,/)  =  0  (2a) 

p(a,l)  =0  (2b) 

the  solution  to  the  Reynolds  equation  (1)  becomes 
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p(x.i)  =  ')x(a-x)  (3) 

where  the  incursion  velocity  v  is  a  function  of  time  i. 

The  normal  force  on  the  blade  tip  due  to  the  pressure 
p(x,r)  is  given  by 

fj=(  (  p(x,t)dxdz  (4) 

Jo  Jo 

Integrating  equation  (4)  with  p(x,t)  given  by  equation  (3) 
yields 

F\=nvb(j\  (5) 

When  the  normal  force  in  equation  (5)  ey  _eds  the  force 
that  can  be  supported  by  the  seal  substrate,  F\  from 
equation  (5)  is  replaced  by  the  following  limiting  elastic 
force  of  the  seal  substrate: 

F\  - t-\i  - C)  (6) 

where  K  is  the  stiffness  of  the  casing. 

The  tangential  force  due  to  shearing  of  the  molten 
metal  is  given  as 

F}=  f  (  T(x,h.r)dxdz  (7) 

J  o  J  0 

where  r(x.h,t)  is  tfie  shear  stress.  The  shear  stress  acting 
on  the  blade  is 


T{xM=T+(V)d£  *iy=h  (8) 

Integrating  relation  (7)  bv  using  equations  (3)  and  (8) 
yields  the  following  ;angenti'>.',  Wade  tip  force: 


(9) 


Abradable  Rub  Model 

Blade  wear,  as  in  the  case  of  the  blade  melting  and 
smearing  on  tl.c  coitiei  substrate,  not  only  produces 
immediate  losses  in  engine  performance  due  to  decreased 
pressure  ratios  but  also  is  a  major  cost  factor  in  engine 
overhaul  since  worn  blades  mun  be  replaced  to  restore 
engine  efficiency.  To  reduce  bl  ide  tip  wear,  abradable 


seal  shrouds  are  being  used.  Current  high-ptessure 
compressor  systems  use  oxidation-resistant  metals  fabri¬ 
cated  to  produce  a  low-strength  abradable  seal,  such  as 
metallic  fibers  or  powders  sintered  into  a  porous,  low- 
density  substrate.  As  the  blade  rubs  on  the  seal,  small 
particles  are  removed  with  each  rub  pass,  analogous  to 
metal  machining. 

Referring  to  figure  2,  the  energy  per  volume  of 
material  removed  per  blade  is  expressed  as 


U  = 


Fi 

(r-C)b 


(10) 


where  r-C  is  the  depth  of  the  cut.  Rearranging  this 
relation  yields 

Fi~  (r-C)Ub  (11) 

A  parameter  study  on  U  was  perfonned. 

The  normal  force  between  the  blade  tip  and  the 
abradable  seal  material  >s  the  same  as  that  given  by 
relation  (6). 


Discussion  of  Results 

The  rotor-bearing  model  described  in  reference  5, 
which  dynamically  simulates  a  typical  small  gas  turbine, 
was  used  as  an  example  problem.  This  rotor-bearing 
model  consisted  of  a  shaft  with  three  disks  mounted  on 
two  axially  stiff  bearings  (fig.  3).  In  this  rotor-bearing 
model  the  bearings  were  mounted  in  squeeze-film  damper 
journals  that  had  centering  springs.  The  rotor  bearing 
model  had  three  undamped  critical  speeds  in  the 
operating  range:  130,  150,  and  190  Hz  (7600,  9200,  and 
1 1  200  rpm). 

The  rotor  was  divided  into  24  segments,  the  same  as  the 
rotor  used  in  reference  5.  The  rotor  was  assumed  to  be 
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'■  Oil  squeeze- 

/  film  damper  Sftai  tub 


balanced  prior  to  the  blade  loss  simulation  and  operating 
at  160  Hz  (9550  rpm).  The  blade  loss  was  simulated  by  an 
instantaneous  application  of  130-^m  (5-mil)  mass 
eccentricity  in  tne  outside  disk  at  the  left  end  of  the  shaft. 
The  blade  tip-,ea!  rub  was  simulated  by  surrounding 
each  disk  with  a  seal  that  had  a  50-^m  (2-mil)  radial 
clearance  and  a  20-MN/m  (100  000-lb/in)  stiffness.  The 
disk  had  rectangular  blades  that  were  13  mm  (500  mils) 
wide  and  2.3  mm  (90  mils)  thick.  The  blades  had  a  tip 
velocity  of  100  m/s  (3800  in/s)  and  only  one  blade  was 
rubbing  at  a  time. 

To  predict  the  rotor  displacement,  the  transient  re¬ 
sponse  code  requires  the  various  interaction  forces.  The 
interaction  between  the  blade  tip  and  sea!  was  calculated 
by  either  the  smearing  or  the  abradable  rub  models 
described  in  the  analysis.  Experimental  data  from  steady 
in-house  rub  tests  seem  to  indicate  for  the  smearing 
model  a  film  thickness  of  2  ^m  (0.08  mil)  and  a  liquid- 
metal  viscosity  of  10  to  50  mN  s/m*  (1.5x10-*  to 

7.5  x  10  - 6  lb  s/in2).  Data  for  the  abradable  model  seem 
to  indicate  an  energy  per  volume  of  material  removed  of 

3.5  to  7.0  GN/m2  (500  000  to  1  000  000  lb/in2).  This 
seemed  to  be  physically  an  order  of  magnitude  too  large, 
since  the  yield  strength  for  most  materals  is  much  less 
then  these  values.  Therefore  this  range  was  reduced  by  a 
factor  of  10.  Four  rub  simulations  were  made.  Two  rubs 
were  simulated  for  each  rub  model.  The  minimum  and 
maximum  of  the  range  of  parameters  given  above  were 
used  to  bracket  the  tesults. 

The  results  of  each  simulation  are  displayed  in  three 
formats.  The  first  is  a  plot  of  the  amplitude  of  the 
displacement  of  the  center  disk  as  a  function  of  time  after 
the  blade  loss.  The  second  is  the  orbit  of  the  center  disk 
after  the  blade  loss.  The  third  is  an  oblique  view  of  the 
envelope  the  centerline  of  the  rotor  swept  out  for  various 
rotations  after  a  blade  loss  had  occurred. 

The  first  model  simulated  was  the  smearing  model.  The 
results  for  the  low-viscosity  case,  p  =  10  mN  s/m2 
(1 .5  x  10-*  lb  sec/in2),  are  shown  in  figures  4  to  6.  The 
results  for  the  high-viscosity  case,  ^  =  50  mN  s/m2 
(7.5  x  10  -*  lb  sec/in2),  are  shown  in  figures  7  to  9, 

In  general,  for  the  low-viscosity  case  the  rotor  seemed 
to  spiral  out  so  that  the  blades  contacted  the  seal  and  then 
began  bouncinR  off  the  seal.  Figure  4  shows  that 


bouncing  ceased  after  25  ms  and  the  blade  tips  w-ere  in 
light  continuous  contact  with  the  seal.  The  orbitai 
response  became  circular  as  shown  in  figure  5.  The  orbit 
stabilized  in  the  forward  whirl  direction  to  a  light  rub 
with  a  maximum  interference  of  5  (0.2  mil).  The 

envelopes  of  the  centerline  of  the  rotor  for  the  first  eight 
rotations  arc  shown  in  figure  6. 

The  amplitude  of  the  bouncing  seemed  to  be  composed 
of  two  frequencies  one  about  70  Hz  (4000  cpm)  and  the 
other  about  550  Hz  (33  000  cpm).  These  frequencies  do 
not  correspond  to  the  original  criticals  of  the  rotor- 
bearing  system  given  in  reference  5.  This  is  not  surpris¬ 
ing,  since  the  original  system  was  modified  with  pseudo¬ 
bearings  at  the  rub  locations.  These  pseudobearings 
stiffened  the  rotor  system  only  during  the  time  the  blade 
tips  were  rubbing  the  seals. 

The  response  for  the  high-viscositv  case  (fig.  7)  started 
out  the  same  as  that  for  the  low-viscosity  case.  The  rotor 
spiraled  out  so  that  the  blades  contacted  the  seal  and  then 
began  bouncing  off  the  seal.  At  this  point  the  two  cases 
began  to  differ.  The  bouncing  (rather  than  decaying) 
seemed  to  grow  in  amplitude  and  the  orbits  became  less 
circular.  The  amplitude  of  the  bouncing  reached  a 
maximum  at  25  ms  (corresponding  to  the  fourth 
rotation).  At  this  point  the  orbit  was  almost  linear  as 
shown  in  figure  8  and  it  marked  a  transition  from 
forward  to  backward  whirl.  As  the  rotor  began  to  whirl 
in  the  backward  direction,  the  amplitude  of  the  bouncing 
began  to  decay.  At  40  ms,  most  of  the  bouncing  ceased 
and  the  blade  lips  were  in  continuous  contact  with  the 
seal.  The  rotor  spiraled  out  to  a  maximum  amplitude  of 
1 30  nm  (5  mil)  at  70  ms,  corresponding  to  a  hard  rub  with 
a  maximum  interference  of  80  /im  (3  mil).  This  seemed  to 
indicate  a  finite  limit  to  the  orbit.  The  rotor  centerline  en¬ 
velopes  for  the  first  eight  rotations  for  the  high-viscosity 
case  are  shown  in  figure  9. 

The  frequency  of  the  blade  tips  bouncing  off  the  Seal 
for  the  high-viscosity  case,  for  times  less  than  40  ms,  was 
about  400  Hz  (24  000  cpm).  This  frequency  was  less  than 
the  frequency  for  the  low-viscosity  case,  550  Hz  (33  000 
cpm).  The  reason  for  this  was  that  in  the  low-viscosity 
case  the  blade  tips  were  in  continuous  contact  with  the 
seal  most  of  the  time,  but  in  the  high-viscosity  case,  the 
blade  tips  bounced  off  the  seal  so  that  contact  between 
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Figure  4  -  Amplitude  of  displacement  of  center  diik  as  a  function  of  time  after  blade  loss  for  vmear'ng  model  (low-visoo«iiv  ca«e). 


Figure  5  -Orbit  of  center  disk  after  blade  loss  for  smearing  model 
(low-viscosity  case).  (Foil  scale  equals  75  ntn  (3  mil).) 


the  blade  tips  and  the  seal  was  very  short.  For  the  high 
viscosity  case  and  times  greater  than  40  ms  the  amount  of 
blade  tip-seal  interference  was  large,  causing  a  large 
normal  force  during  the  time  of  rubbing  and  therefore 
generating  high-frequency  vibration  components,  1000 
Hz  (60  000  cpni). 

The  second  model  simulated  was  the  abradable  model. 
The  results  for  the  low-energy-per-volume  case,  (7=350 
MN/m2  (50  000  lb/in2),  are  shown  in  figures  10  to  12. 
The  results  for  the  high-energy-per-volume  case,  (7  =  700 
MN/m2  (100  000  lb/in2),  are  shown  in  figures  13  to  15. 
The  abradable  model  necessarily  started  out  the  same  as 
the  smearing  model.  The  rotor  spiraled  out  so  that  the 
blades  contacted  the  seal  and  then  began  bouncing  off  the 
seal.  At  this  point  the  two  models  began  to  differ. 

The  results  shown  in  figure  10  for  the  low-energy-per- 
unit-volume  case  in  the  abradable  model  are  similar  to 
those  for  the  low-viscosity  case  in  the  smearing  model. 
Comparing  the  orbital  response  for  the  smearing  model 
low-viscosity  case  (fig.  5)  with  the  orbital  response  for  the 
abradable  model  (fig.  11)  shows  that  the  envelopes  that 
the  rotor  centei lines  swept  out  for  the  two  cases  were 
virtually  indistinguishable.  Like  the  low-viscosity  case,  in 
the  abradable  model  the  blade  tips  were  in  light 
continuous  contact  with  the  seal  and  the  orbit  stabilized 
in  the  forward  whirl  direction  to  the  same  light  rub  a 
maximum  interference  of  5  (0.2  mil).  The  amplitude 
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Figure  6  -  Envelopes  of  roior  centerline  for  first  eight  rotations— 
smearing  model  (lo» -viscosity  case).  (Full  scale  equals  130  urn 
(5  mils)  ) 

of  the  bouncing  seemed  to  be  composed  of  the  same  two 
frequencies,  70  and  550  Hz  (7000  and  33  000  cpm). 
However,  the  time  of  decay  was  much  (aster  (30  ms)  for 
the  low-energy  case  than  for  the  low-viscosity  case. 

For  times  less  than  40  ms  the  results  for  the  high- 
energy-per-unit-volume  case  in  the  abradable  model  (fig. 


Figure  8  -Orbit  of  center  disk  after  blade  loss  for  smearing  model 
(high-viscosity  case).  (Full  scale  equals  150  (6  mil).) 


13)  are  similar  to  those  for  the  high-viscosity  case  in  the 
smearing  model  (fig.  7).  The  envelopes  the  rotor  center- 
line  swept  out  for  the  smearing  model  were  similar  to 
those  for  the  abradable  model  with  the  exception  that  the 
normal  force  for  the  smearing  model  was  much  higher. 
The  amplitude  of  the  bouncing  reached  a  maximum  at 
the  same  25  ms,  corresponding  to  the  fourth  rotation.  As 
in  the  high-viscositv  case  this  was  also  the  transition  from 
forward  to  backward  whirl.  The  frequency  of  the 
bouncing  was  the  same,  400  Hz  (24  000  cpm).  This 
indicated  that  the  time  of  contact  between  the  blade  tips 
and  the  seal  was  short.  At  about  40  ms  the  bouncing 
decayed  and  the  blade  tips  were  in  continuous  contact 


Time,  ms 

Figure  7.  -  Amplitude  of  displacement  of  center  disk  as  function  of  time  after  blade  loss  for  smearing  model  (high-viscosity  case). 
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gurc  12  -  tmcliipei  of  rolor  ccmcrlinc  lor  tirsi  ciglil  rolations — 
abiadahic  model  (lo»-eucrg>  easel.  (I  ull  stale  equals  1)0  <on 
()  mill  > 


with  the  seal.  At  this  point  the  two  models  began  to 
differ.  The  rotor  began  to  spiral  out  and  the  amplitude 
seemed  to  grow  without  limit,  as  shown  m  figures  1.)  and 
14.  The  rale  of  growth  was  very  fast.  This  was  similar  to 
the  dry  friction  model  (ref.  5)  when  the  coefficient  ol 


f- 


a — j- 


Figure  M  •  Orbn  of  center  disk  alter  blade  loss  for  abradable  model 
(litgh-ett;*rg>  cave)  (f  ull  scale  equals  600  ,«n»  (24  mil).) 


y 


figure  15  1-mdopcx  of  rotoi  centerline  lor  fir  si  sesen  roiations— 
abradable  model  (high-energx  case)  (lull  scale  equals  130  >*m 
(.<  mil).) 


friction  exceeded  a  particular  threshold.  The  rotor 
centerline  envelopes  for  the  loss-  and  high-energy 
abradable  models  are  shown  in  figures  12  and  15, 
respectively. 

In  general,  by  adjusting  the  patameters  used  m  the 
model,  it  was  possible  to  produce  similar  results  with  any 
model  up  to  the  point  at  which  the  rotor  went  into 
backward  whirl.  Once  the  rotor  went  into  backward 
whirl,  the  type  of  model  used  was  important. 


Conclusions 

A  study  was  conducted  of  the  response  of  turbine 
engine  rotors  to  two  different  blade  tip-seal  interference 


rub  models.  The  following  conclusions  were  drawn: 

1.  The  abradable  model  was  more  sensitive  to  small 
changes  in  the  energy  per  unit  volume  ol  material 
removed  than  the  smearing  model  was  to  changes  in  the 

molten  metal  viscosity. 

2.  Both  the  abradable  and  smearing  models  had  a 
threshold  which  when  exceeded  caused  the  rotor  to 
proceed  into  backwaid  whirl. 

3.  When  the  abradable  model  went  into  backward 
whirl,  n  resulted  in  a  catastrophic  failure  whereas  the 
smearing  model  resulted  in  a  more  benign  failure. 

In  general  the  two  models  can  be  manipulated  to 
produce  similar  results  up  to  the  time  at  which  the  rotor 
goes  into  backward  whirl.  After  this  time  the  type  of 
model  is  important. 

National  Aeronautics  and  Space  Administration 

Lewis  Research  Center 

Cleveland,  Ohio,  September  27,  1973 
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ie  Abstract 

Using  a  direct-integration,  transient-response  rotor  dynamics  computer  code,  the 
response  of  turbine  engine  rotors  to  two  different  blade  tip  -  seal  interference 
rub  models  was  studied.  The  first  model,  an  abradable  seal  rub  model,  is  based  on 
an  energy-loss-per-unit-volume  theory.  It.  is  applicable  to  a  ceramic  turbine  blade 
tip  seal.  The  second,  a  smearing  model,  is  based  on  viccous  hydrodynamic  theory. 

It  is  applicable  to  a  metallic  blade  tip  seal.  The  results  from  these  two  models 
were  compared  with  those  from  a  previously  studied  model  based  on  dry  friction 
theory.  The  abradable  model  was  very  sensitive  to  small  changes  in  the  energy  per 
unit  volume,  and  once  a  threshold  was  exceeded,  the  rotor  went  into  a  backward 
whirl.  The  amplitude  seemed  to  grow  without  limit.  This  was  similar  to  the  dry 
friction  model  when  the  coefficient  of  friction  exceeded  a  particular  threshold. 

The  smearing  model  was  not  as  sensitive  to  small  changes  in  the  viscosity,  but  a 
threshold  viscosity  was  found.  When  it  was  exceeded,  the  rotor  went  into  backward 
whirl,  but  the  amplitude  seemed  to  grow  to  a  finite  limit.  - 
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